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Abstract 

The  covariance  function  of  the  Newtonian  potential  of  a  random  orthogonal 
I 

set  function  on  the  unit  sphere  in  three  dimensions  is  derived,  and  it  is  shown 

that  the  coefficients  to  the  series  expansion  of  this  are  simply  related  to  the 

moments  of  the  covariance  measure  of  the  random  set  function. 

Furthermore,  as  an  application,  it  is  shown  that  available  gravity  data 

indicates  a  mass  distribution  inside  the  Earth  which  becomes  more  and  more 

I 

irregular  as  one  approaches  the  center  of  the  Earth. 
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1.  Introduction 


In  Lauritzen  (1973)  a  derivation  of  stochastic  models  for  the  gravity 
potential  of  the  Earth  was  given.  It  was  shown  how  measurements  of  gravity 
anomalies  and  deflections  of  the  vertical  can  be  treated  as  linear  operations 
on  the  random  potential  function,  thus  giving  opportunity  to  estimate  the 
covariance  function  under  the  assumption  that  it  be  invariant  under  the  ortho¬ 
gonal  group,  and  to  predict  linear  functionals  of  the  potential  where  these 
are  unknown.  A  specific  expression  of  the  covariance  function  was  shown  to 
fit  data  very  well. 

However,  the  derivation  was  based  entirely  on  analytical  properties  of 
harmonic  functions  and  did  not  utilize  the  relation  between  the  mass  density 
in  the  interior  of  the  Earth  and  the  corresponding  potential,  thus  making  it 
difficult  to  give  any  physical  interpretation  of  the  results. 

The  present  paper  gives  a  derivation  of  stochastic  models  for  the  dis¬ 
turbing  potential  based  on  the  variations  of  the  mass  and  it  is  shown  that 
the  available  data  indicates  the  mass  distribution  becoming  more  and  more 
irregular  as  one  approaches  the  center  of  the  Earth. 

The  disturbing  mass  is  described  by  a  random  orthogonal  set  function 
and  it  is  shown  that  the  "degree-variances"  of  the  potential  covariance 
function  are  simply  related  to  the  moments  of  the  covariance  measure  of  the 
random  orthogonal  set  function. 

As  a  consequence,  it  is  shown  that  the  problem  of  determining  the  possible 
covariance  functions  for  the  disturbing  potential  arising  from  models  with  the 
above  mentioned  structure  is  equivalent  to  the  classical  Hausdorff  moment  problem. 


Finally,  the  covariance  measure  corresponding  to  the  model  proposed 
in  Lauritzen  (1973)  is  derived. 


2.  Potentials  of  random  orthogonal  set  functions  on  the  unit  sphere 
in  3-dimensional  Euclidean  space. 

Let  Z  be  a  random  orthogonal  set  function  on  the  Borel  sets  of  S  , 

3 

the  closed  unit  sphere  in  R  ,  3-dimensional  Euclidean  space,  with  covariance 
measure  y  satisfying  y(s)  <  +oo  ,  i.e.  Z  is  a  stochastic  process  on  the 
Borel  sets  of  S  with 

i)  EZ (A)  =  0 

ii)  EZ (A)Z (B)  =  y(A  0  B) 

iii )  Z  (A)  +  Z  (B )  =  Z  (A  U  B)  a.s.  for  A^l  B  =  0  ,  A  and  B  being 
arbitrary  Borel  sets  of  S  . 

For  t  $  S  ,  define  the  stochastic  integral  in  the  L^-sense  by 


£(t)  = 


u6S 


t-u 


Z  (du) 


This  defines  a  stochastic  process  on  R^VS  with  mean  value 


E£(t)  = 


u£S 


,  — "I'T  EZ  (du)  =  0 
t-u 


and  covariance  function 


E5(s)£(t)  = 


f 


u^S 


T-n: — r  T- — rr  E{z(du)z(dv)} 
v£S  TFFTT|s-v 


u£S 


t-u  s-v 


y (du)  , 


see  e.g.  Grenander  and  Rosenblatt  (1957)  p.  25  ff. 


2 


The  stochastic  process  E,  shall  be  called  the  Newtonian  potential  of 

Z  ,  in  accordance  with  the  classical  potential  theory,  where  the  Newtonian 

3 

potential  of  a  Radon  measure  v  on  a  compact  set  KC  R  is  defined  as  the 

"5 

function  4> :  R  \K  -*■  R  given  by 


(p  U) 


r 


y6K 


As  we  have  the  expansion  (e.g.  Hobson,  1955) 


where  ^  is  the  angle  between  the  t  and  the  u  vectors  and  Pn  are 
the  Legendre  polynomials,  we  get 

f  lluMn+lnpn(COS  *+JPm(cOS 


E£(s)£(t)  = 


oo  oo 

I  I 

n=0  m=0 


ues 


ihir 


I® 


y(du)  . 


If  we  assume  that  y  is  uniform  on  spheres  of  the  form  {| |u| |  =  p},  and 
introduce  spherical  coordinates 


s  = 


/  sl\ 
S2 
S3  / 


/r s  sin  0S  cos  Xs\ 
=  I  r  sin  0  sin  X 


l*s  cos  9s 


rM 

— 

6s 

L  X  / 

/ 

\  s  / 

we  get 

E5(s)g(t) 


OO  oo 

-i-  y  i 

rs  rt  m=0  n=0\ 


1  n+m 
x 

.  n  m 
Or  r, 
s  t 


f 

2lT 

l  i 

0 

/  VJ 

* 

a  (s,t,0,A)sin  0d0dX 
nm  ’  ' 


where 
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y(dr,d0,dA)  =  ri(dr)  sin0  d0  dA 


and 


a  (s,t,0,A)  =  P  (cos  )  P  (cos  \p,  )  ,  u  ^  (l  0  A)  . 
nm  n  su  m  tu  ’  ’ 


Using  the  well-known  formula 


Pn(c°s  ipsu)  =  Pn(cos  0s)  Pn(cos  0)  +  2  l Js_0|pni(c0s  0s)  Pni(cos6)cos  i^-A) 


i=l 


where  P  ^  are  the  Legendre  functions,  we  get 


nm 


n  m 

(s,t,e,A)  =  l  l  b  (s,t,0,A) 
i=0  j=0  J 


where 


b.  (s ,t ,0  ,A )  =  c.  .  [R  . (0  ,A  )R  .  (0.  ,A,  )R  . (0,A)R  .(0,A) 
ljnm  ’  ’  ’  '  ljnm  ni  s’  s  mj  t’  t  m  *  mj  ’ 


Here 


+  R  .(0  ,A  )S  . (0 .  , A  )R  .(0,A)S  .  (0,A) 
m'  s’  s  mj  t ’  t '  ni  ’  mj 


+  S  .  (0  ,A  )R  .  (0  A  )S  .  (0 ,  A  )R  .  (0  ,A ) 
ni  s  ’  s  mj  t  ’  t  ni  mj 


+  S  . (0  ,A  )S  . (6  ,A  )S  . (6,A)S  .(6,A)] 
ni  s’  s  mj  t  ’  t  ni  ’  mj  ’ 


R  .  (Q,A)  =  P  .  (cos  0)  cos  iA 
ni  ni 


S  .  (0,A)  =  P  .  (cos  0)  sin  iA 

ni  ni 


are  the  surface  spherical  harmonics.  Using  the  orthogonality  property  of 


k 


these,  the  only  terms  in  the  sums  contributing  anything  to  the  integral  will 


he  those  of  the  form  R  .R  .R  .R  .  or  S  .S  .S  .S  .  .  As  for  i  4  0 

ni  ni  ni  m  ni  ni  ni  ni 


'27T 

o 

’2TT 

0 

R  .  (0,A)  sin9  d0  dA  = 

Ill 

0 

0 

7T 


52.  (9, A)  sin6  d9dA  =  , 

ni  5  (2n+l ) (n-i ) !  5 


f  2tt 
0 


rir 


R2o(0,A)  Sin6  d9  aA  =  2S1  ’  Sni  =  0  ’ 


and 


f  l  J  (n~i  )  •  \  2 


Trv+ i  ) ! 


C.  . 

linn 


J 


for  i  4  0 
otherwise , 


we  obtain  by  using  (*)  again,  that 
f 27T  fTT 


0 


0 


a  (s,t,6,X)  sin0  d9  dA  =  6  - 

nm  nm  2n+l  n 


P  (cos  ) 
”  st7 


Hence 


,  1  ,ntl  Pn(cos  *Sjt) 


I  (rV) 

n=0 \  s  t / 


x2nri  ( dx ) 


Defining  an  = 


x2  n  nM  ,  we  obtain  that 


a  =  a  *  y 
n  rn 


where  v  is  the  n’th  moment  of  the  distribution  of  the  random  variable  X  , 
'n  ’ 


where  X  follows  the  distribution 


rx 


P{X  <  x}  = 


n(<iu) 


n(du) 
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In  terms  of  the  covariance  measure ,  we  get 


u  y  (du)  . 


The  problem  of  when  a  prescribed  sequence  can  be  a  sequence  of  moments 

of  a  distribution  on  the  unit  interval  is  the  classical  Hausdorff  moment  problem, 

having  the  solution  that  (Y  ,  n=0,l,2,  ...  )  is  a  moment  sequence  iff  it  is 

completely  monotone,  i.e.  iff  (-l)VAVy  >  0  for  all  v,n  ,  where  Ay  =  y  _  y 

'n  =  n  n+1  n 

see  Feller  (1966  ) . 


3.  Application  to  the  gravity  potential  of  the  Earth. 

Lauritzen  (1973)  studied  stochastic  models  for  the  disturbing  gravity 
potential  of  the  Earth,  and  showed  that  the  general  form  of  a  rotational  in¬ 
variant  covariance  function  for  a  random  potential  is 

00  p  /  1  \n+1 

s)£(t)  =  I  an  -  Pn  (cos  0)  , 

n=0  \  s  t / 


where  a^  just  have  to  be  so  that  the  series  is  convergent. 

The  random  disturbing  potential  can  be  represented  in  a  series  expansion 
with  random  coefficients  a  and  8 


5(t)  = 


OO  ••  T 

i  fH 

n=0  At  ) 


nm 
■v  n+1  /  n 


m=0 


nm 


cx  R  _  (0,  ,X.  ) 
nm  nm  t  t 


n 

Y  8  s  (e.,x.) 

L  nm  nm  t  t 
m=0 


As  one  has  rather  precise  knowledge  of  the  coefficients  ^nm5  3nn  for  n=0,l,2, 
one  could  consider  those  as  fixed,  i.e.  consider  the  covariance  function  of 
E{ t)  conditional  on  the  known  values  of  a  ,  6  for  n  =  0,  1,  2,  as  also 
done  in  Lauritzen  (1973);  this  can  be  shown  (same  paper)  to  affect  the  covariance 
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2  2  2 

function  only  by  putting  =  a2  =  0  . 

From  an  empirical  covariance  function  obtained  from  gravity  anomalies, 
it  was  demonstrated  that  an  expression  of  the  form 


for  n  =  0,  1,  2 


forni3 

gave  an  extremely  good  fit  to  the  observations. 

The  investigations  in  section  2  of  the  present  paper  give  us  a  possibility 
of  an  interpretation  of  these  coefficients  in  terms  of  the  variations  in  the 
mass  distribution  inside  the  Earth.  We  have  the  relation 


2  2  kir  2  2n+l  2 

n  n  2n+l  n  4  if  n 

Now,  suppose  that  the  Earth  is  the  unit  sphere  in  F3  and  the  deviations 
between  a  homogeneous  and  the  actual  mass  distribution  is  given  by  a  random 
orthogonal  set  function  Z  with  covariance  measure  absolutely  continuous  w.r.t. 
Lebesgue  measure  on  S  and  density  f  given  by 


f(u) 


for  | |u| |  <  e 


e  <_  |  |  u|  |  <^1  . 


This  corresponds  to  p(dx)  =  g(x)dx  ,  where 


for  x  <  e 


for  e  <_  x  <_  1  , 
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corresponding  to  for  n  >_  3 


2  _  _A  2n+l  _A  2n-U  /  3ef_  _5_  \ 

°n  "  4tt  (n-1 )  (n-2 )  Utt  £  ^n-1  '  n-2  J 


For  £  being  small  we  then  have  the  approximation 


2  A 

an  '  (n-1)  (n-2) 


+  0(£2)  •  O(i-)  for  n  >  3  , 


the  main  term  of  which  exactly  is  the  one  proposed  by  Lauritzen  (1973). 

2 

The  good  fit  of  this  form  for  a^  should  then  indicate  a  mass  dis¬ 
tribution  inside  the  Earth  which  gets  more  and  more  irregular  in  absolute 
variation  as  one  approaches  the  center  of  the  Earth. 
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